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1. Introduction 

Energies and transition probabilities of K/5 hypersatellite lines are evaluated in this work 
for selected atoms throughout the periodic table. In a previous paper [T|, we reported 
on calculated values of Ka hypersatellite line energies for atoms with 12 < Z < 30. 

A hypersatellite line is an X-ray line for which the initial state has two vacancies 
in the same shell. This is the case, for example, when a double ionized K-shell state 
decays through the transition of one M-shell electron. Lines corresponding to ls~^ —>■ 
ls^^3p^^ transitions, where ni~^ means k vacancies in the ni subshell, are called K/Jj^g 
hypersatellite lines. 

Atoms where a whole inner shell is empty while the outer shells are occupied were 
first named hollow atoms by Briand et al |2j and are of great importance for studies 
of ultrafast dynamics in atoms far from equilibrium and with possible wide-ranging 
applications in physics, chemistry, biology, and materials science j^]. 

Very scarce experimental data exist on energies of K hypersatellite lines, due to 
the low probability of creation of the initial two K holes. Briand et al ^ used a 
coincidence method to study the K hypersateUite spectrum in Ga and measured a 390 
eV energy shift for the K/?*^ lines relative to the corresponding diagram lines. Energies 
of K/3^ hypersatellite lines were later measured by Briand et al |2l El for Mn, Ga and 
In. Diamant et al jZj obtained high resolution K/?|^3 and Ka*^ hypersatellite spectra 
of Fe, using monochromatized synchrotron radiation photoexcitation. Similar work has 
been performed for Cr K/5^ hypersatellite spectra jH] . More recently, this work has been 
extended to Ti 0. 

On the theoretical side, calculations of the energies and transition probabilities 
of K/5j^3 hypersatellite lines have been performed [10] using a Dirac-Hartree-Slater 
approach. This approach employs a local approximation to the atomic potential. With 
the wave functions obtained from this potential, perturbation calculations were then 
used to obtain the energies. Diamant et al |7] performed relativistic multi-configuration 
Dirac-Fock calculations for the ¥^0^^ hypersatellite lines of Fe to compare with their own 
experimental findings. 

The energy shifts between the K/^j^g hypersatellite lines and the corresponding K/?i 3 
diagram lines (as the Ka^-Ka ones) are considered to be very sensitive to the Breit 
interaction. Chen et al JU] calculated the influence of the Breit interaction on these 
shifts, using a mathematical form appropriate for the local approximation employed. 
To check the importance of the Breit interaction, values of these shifts as well of the 
K/Sj^/K/^g intensity ratios were calculated in this work for selected values of Z, with and 
without inclusion of the this interaction. 

Using the results of our calculations, we were able to obtain the widths of the ls~^ 
two-hole levels of Al and Sc. 
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2. Calculation of atomic wave functions and energies 

Bound state wave functions and radiative transition probabilities were calculated 
using the multi-configuration Dirac-Fock program of J. P. Desclaux and P. Indelicato 
im IT2| IT3J . The program was used in single-configuration mode because correlation 
was found to be unimportant for the energies and probabilities transition. The wave 
functions of the initial and final states were computed independently, that is, atomic 
orbitals were fully relaxed in the calculation of the wave function for each state, and 
non-orthogonality was taken in account in transition probability calculations. 

In order to obtain a correct relationship between many-body methods and quantum 
electrodynamics (QED) JEIEIQSIEI; one should start from the no-pair Hamiltonian 

N 
i=l i<j 

where r^j = |rj — Tj|, and Tin is the one electron Dirac operator and V is an operator 
representing the electron-electron interaction of order one in a, properly set up between 
projection operators yl^^ = ^f ■^'j ^^ avoid coupling positive and negative energy states 

V(r.,) = Al^Vir^^Al^. (2) 

In nonrelativistic quantum mechanics the interaction between charged particles is 
described by an operator associated to the potential energy of interaction, V {rij), 
which is only a function of the distance between particles, Tjj. However, in quantum 
electrodynamics there exists no potential energy of interaction that depends on the 
coordinates of the interacting charged particles taken at the same time, since the charges 
interact not directly, but via the electromagnetic field. The operator for the interaction 
potential energy is replaced, in the first approximation, by the scattering matrix |TH] 

^%^cD = ^'^ H^D (y) 7mV^b (y) D^ (y - ^) V^c (y) 7/.^a (y) d'^y d'^x, (3) 

After performing the time integration of (jH)), one finds the following well known 
interaction potential, in Coulomb gauge and in atomic units, 

9 (rij) = — (4a) 



OLi ■ Otj 



m 



COS ^^^^^ - 1] 



fij ^ C 

COS (iiia^) - 1 
+ K ■ V,)(a, ■ V,) ^-^ , (4c) 

where is Uij is the energy of the photon exchanged between the two electrons, ckj are 
the Dirac matrices and c = 1/a is the speed of light. The first term ()4aj) represents 
the Coulomb interaction, the second one (J4S|) is the Gaunt interaction (magnetic 
interaction), and the last two terms (|¥c|) stand for the retardation operator. In this 
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expression the V operators act only on Tjj and not on the following wave functions. By 
a series expansion of the terms Eqs. (J4^-(0cj) in powers of Uijrij/c <^ 1 one obtains the 
Breit interaction, which includes the leading retardation contribution of order 1/c^. The 
Breit interaction is the sum of the Gaunt interaction (J4^ and of the Breit retardation 



9^ in,) = ^^^ - -- -^^ -^^ (5) 

In the present calculation the electron-electron interaction is described by the sum of 
the Coulomb and the Breit interaction. Higher orders in 1/c, coming from the difference 
between Eqs. Pcj) and (0) are treated here only as a first order perturbation. 

We use the Coulomb gauge as it has been demonstrated that it provides energies 
free from spurious contributions at the ladder approximation level and must be used in 
many-body atomic structure calculations J19j . 

Finally, from a full QED treatment, one also obtains the radiative corrections 
(important for the innermost shells) to the electron-nucleus interaction (self-energy and 
vacuum polarization). The one-electron self-energy is evaluated using the one-electron 
values of Mohr and coworkers |2II1 El 122] • The self-energy screening is treated with the 
Welton method developed in Refs. j2Sl 1211 123 I2HI • This method yields results in close 
agreement (better than 5%) with ab initio methods based on QED j2Zl I2H1 12^1 1^ - 
without the huge amount of effort involved. 

The vacuum polarization is evaluated as described in Ref. [HI]. The Uelhing 
contribution is evaluated to all orders by being included in the self-consistent field 
(SCF). The Wichmann and Kroll, and Kallen and Sabry contributions are included 
perturbatively. These three contributions are evaluated using the numerical procedure 
from Refs. JS21ISS]. 

3. Results 

3.1. Introduction 

We calculated the energy shifts of the K/^j^g hypersatellite lines relative to the 
corresponding diagram lines for several atoms with 13 < Z < 80, and the Kj3i and 
K/^g energy shifts, and their relative intensities, for selected atoms with 18 < Z < 80. 
The wave functions of the initial and final states were computed independently, that is, 
atomic orbitals were fully relaxed in the calculation of the wave function for each state. 

In what concerns the precise identification of the K/3j^3 hypersatellite lines some 
comments are in order. 

Depending on the configurations of the initial and final states, for the different 
values of Z, the number of transition lines that must be computed may range from only 
two, when the initial state has only closed shells, to several thousand, when unfilled 
shells exist. 

For elements with filled shells, namely Ar, Ca, Zn, Kr, Sr, Pd, Cd, Xe, Ba and Hg 
the ls~^ ground configuration corresponds to only one level, the ^Sq level, and each of 
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the Kf]^ lines is identified by a precise level transition, 
K/Jg^: ls-2 iSo -^ ls-i3p-i iPi 
K/?^ ls-2 iSo ^ ls-i3p-i 3p^. 

By analogy with the corresponding diagram lines, the line that corresponds to the 
larger transition energy value is labeled Kf3^ and the other one is labeled Kj3^. 

For Z > 29 there are two sets of transition lines, separated in energy by more than 
3 eV, corresponding to the K(3^ and K/^g lines. In the first set the decay is due mainly 
to the 3p3/2 electron transition, whereas in the second set is mainly due to the 3pi/2 
electron transition. LS coupling dominates the level structure for elements with Z < 29, 
as can be seen in Fig. ^for Si. In this figure, the labels refer to the only transitions that 
contribute for the spectrum. Intercombination lines give negligible contribution. 

3.2. Transition probabilities 

The transition probability W^ for the line X is defined as 

where 7 is, in this case, a given double-K hole configuration, A^ (i) is a collection of 
excited atoms per unit volume and N (7) = ^j A^(i). W^ is the transition probability 
for the line X from an initial level i, defined by 

W,'' = J2w.f, (7) 

/^ 

where f^ runs over all possible final levels in the radiative de-excitation process leading 
to the X line, and Wij is the probability per unit time that an atom in the excited level 
i, will decay to the level /, through the spontaneous emission of a photon. 

For short lifetimes r of the excited levels, compared with characteristic creating 
times (the inverse of the number of excitations per second undergone by one atom), 
the number of atoms doubly ionized in the K shell created in the excited level i per 
unit time C^ (i) equals the rate at which the atoms leave the level i, N (i) /r (i), by all 
possible transitions. Assuming that the i level of a given double-K hole 7 configuration 
is fed according their statistical weight, we have 

C,(z) = C,^, (8) 

where g (i) and g (7) are the multiplicities of the i level and of the 7 double-K hole 
configuration, respectively, and C^ is the number of double-K ionised atoms created per 
unit time and per unit volume. From Eq. (jH)) and (jH)) we obtain 



rx J:^9{^r{^)W,^ 



Using 



w^ = ^'''\.:\' ■ (9) 



'^^ Z,w,.Zr^,r (^°' 
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where / runs for all possible final levels that can be reached by radiative transitions, 
with probabilities Wif, and /' runs for all possible final levels that can be reached by 
radiationless transitions, with probabilities Aifr, we get 

W^ = ^'+^' ■ (11) 



E. 



5(0 



We made use of Wi = J2f ^if^ ^^^ Ai = J^f Aif. 

For the elements where the two K-hole ground configuration has more than one 
level, it is therefore necessary to compute, for each of those levels, not only radiative 
transition probabilities, but also the radiationless (Auger) transition probabilities, to 
obtain the quantities W^ . 

A complete calculation of radiative and radiationless decay rates from the double-K 
hole ground configuration was performed for Al and Sc. Radiative transitions include 
Ka^ and Kf]^ hypersatellite lines, as well as Kaa (one electron - two photon transitions) 
lines. The results are presented in Table H together with the total radiative and 
radiationless transition probabilities, for each of the two initial levels. 

In Table |2we provide the results of the complete calculation of W^ for the different 
lines in Al and Sc and of the statistical average transition probability of line X, defined 
as the quantity 

^sA = ^^J:9i^)W- (12) 

We observe that the values of W^ and Wgj^ are nearly identical. This results from 
the fact that Eq. ()12j) can be obtained from Eq. (fTTjl if the summation Wi + Ai has the 
same value for all initial levels. This is the case of Al and Sc presented in Table ^ 

Total radiationless level widths are the sums of a large number of transition rates. 
We may assume that the relativistic effects tend to average out to some extent, similarly 
to what happens with total radiative widths |S1]. To test this assumption, the value 
of radiationless (Auger) decay rate for each initial level of the Ti 2s^ 2p^ 3s^ 3p^ 3d^ 
4s^ ground configuration was computed by adding the values of radiationless transition 
probabilities for all levels of the final Is 2p^ 3s^ 3p^ 3d^ 4s^ configuration. As shown 
in Table El no significant variation of the decay rates was found for different initial 
levels, which shows that the radiationless decay rates do not depend significantly on the 
particular level of the initial configuration. This validates the use of Eq. (II 2j) as a good 
value to W^. 

The ratio of the intensities of K/3j^ to K/Jg hypersatellite lines computed in this 
work, with and without the Breit interaction, using the MCDF code of Desclaux and 
Indelicato, are presented in Table IH together with the values obtained by Chen et al 
, which included Breit and vacuum-polarization corrections. We notice that our values 
for these ratios are larger than Chen's results for Z < 40. The two approaches yield 
ratios that are in good agreement for Z > 40. 
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In order to compare the transition energy values obtained in this work with 
experiment and calculations from other authors, we used the statistical average energy 
-EsA fo^ ^he X line defined in our previous article pj as 

(13) 



^SA= ^-T^^W 



/ 



X 






In this calculation we assumed that all i levels of the 7 configuration are statistically 
populated. The quantity in parenthesis is the average energy of the X line, defined as 
the sum of the energies of all individual i ^ f transitions in the X line from an initial 
level i, Eif, weighted by the corresponding Wif radiative transition probability. 

3.3. Widths of Is'"^ levels of Al and Sc 

Using the values presented in Table QJ we were able to calculate the widths of the ls~^ 
levels of Al and Sc, which are displayed in Table El 

We beheve this is the first time that level widths are calculated for a double-K hole 
level. These values can be compared with existing single-K hole level widths, using the 
expression T (ls~^) = 2T (ls~^). 

For Al, using the Evaluated Atomic Data Library (EADL) value T (is^^) = 0.37 
eV proposed by Campbell and Papp [33], we obtain T (is^^) = 0.74 eV, lower than 
the value calculated in this work. We note, however, that the experimental values of 
r (ls~^) for Al referred to by the same authors are higher than their proposed value. 
Three of these values were derived from indirect measurements, using other level widths 
to obtain the ls~^ level width. The only experiment that led directly to the ls~^ level 
width yielded the value 0.47 eV. 

For Sc, the same authors propose F (is^^) = 0.83 eV, which yields F (is^^) = 1.66 
eV, in excellent agreement with the value obtained in this work. 

3.4. Energy shifts 

In Table El we present the results obtained in this work for the Kj3i and K/Jg energy 
shifts of the elements where we can distinguish these two lines. 

This Table shows that our results for the KjS^ hypersatellites energy shifts, relative 
to the corresponding diagram line energies, are in good agreement with the results of 
Chen et al \TU\, ours being smaller by less than 0.2 % throughout. 

To compare to the available experimental results, we present in Table [7|the K/5|^3 
energy shifts calculated in this work. Our results agree in general with experiment, as 
it can be seen in Fig. although the uncertainties of the latter are very large, with the 
exception of the recent experimental value of Diamant et al ^ . 
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3.5. Breit interaction and QED corrections 

To assess the contribution of the Breit interaction to the K/?'^ hypersatelhtes energy 
shifts, we computed these shifts with and without inclusion of the Breit interaction 
in the calculation. We computed separately the K/3|^ energy shifts, first with the Breit 
term (cf. Eq. ^) included in the self-consistent process and the higher-order terms as a 
perturbation after the self-consistent process is finished, and then with Breit interaction 
neglected. The results are presented in Table |H1 Although Chen et al [TU1 present 
their results for these shifts, obtained with the Dirac-Slater approach, in graphic form 
only, we easily conclude that our results, using the MCDF approach, are in very good 
agreement with the results of Chen et al . 

The K/5j^ to K/33 intensity ratio is sensitive to the inclusion of the Breit interaction, 
similarly to Chem et al ^0] finding for the ¥^a\ to Ka| intensity ratio. The inclusion 
of this interaction decreases the 'KP^ intensity ratio at low Z (21% for Z = 18) and 
increases it for medium and high Z (~ 5% at Z ~ 50). However, since relativity affects 
the 3pi/2 and 3p3/2 levels in a similar way the KjS^ intensity ratio increases monotonically 
towards the jj coupling limit of 2. 

The evolution of the relative contribution of Breit interaction to the Kf3^/Kf3^ 
relative intensity ratio is illustrated in Fig. El 

On the other hand, QED contributions for the energy shifts and transition 
probabilities have been found to be negligible. For instance, QED contributions for 
the K/3j^ — K/5i energy shift in Hg is only 0.3% compared with 13% from the Breit 
interaction contribution, whereas for the Kf3^/K(3^ intensity ratio the QED contribution 
is 0.05%, compared with 1.8% from the Breit contribution. The QED contributions for 
the Kj3i — K/3i energy shift are presented in Table |H1 

4. Conclusion 

In this work we used the MCDF program of Desclaux and Indelicato to compute energy 
shifts of the K(3^ and K/^g hypersatellite lines relative to the parent diagram lines for 
several values of Z throughout the periodic table. One of the aims of this work was to 
assess the contribution of the Breit interaction to these shifts. Our results confirm the 
earlier findings of Chen et al fOI for these shifts and extended them to higher values 
of Z. We also calculated the Kj3i to Kj3^ intensity ratio for the same values of Z. Our 
results are significantly lower than Chen et al values for the same ratios, for Z lower 
than 40, and agree with the values of these authors for higher values of Z. The total 
widths of the double-hole K levels of Al and Sc were also computed and our values were 
found in good agreement with the ones obtained from proposed values of single-K hole 
levels 
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Table 1. Radiative transition probabilities of Kaa, Ka*^ and K/?*^ lines, Auger and 
total transition probabilities for each initial level (LSJi) in Al (Z = 13) and Sc 
(Z = 21). Numbers in parenthesis indicate a power of ten. 









Wi" (s-^) 






W, (s-i) 


A. (s-^) 


LSJ, 


Ka^ 


Ka^" 


K/3S>3 


Ka2a3 


Kaias 




Al 2p^/2 

'P3/2 
Sc 2D3/2 


7.35(13) 
6.36(13) 
4.93(14) 
4.88(14) 


6.74(11) 
5.40(11) 
2.98(13) 
2.18(13) 


1.20(12) 
1.20(12) 
7.04(13) 
7.01(13) 


9.37(10) 
1.00(11) 
4.01(11) 
3.92(11) 


5.45(7) 
5.06(7) 
1.69(9) 
1.19(9) 


7.55(13) 
6.55(13) 
5.94(14) 
5.81(14) 


1.56(15) 
1.49(15) 
1.91(15) 
1.95(15) 



Table 2. Comparison between the results of a complete calculation W^ and a 
statistical average calculation VFgA of the transition probability for Kaa, Ka^ and 
K(3^ lines in Al {Z = 13) and Sc {Z = 21). Numbers in parenthesis indicate a power 
of ten. 

Ka^ Ka^l K/Jj'a Kaa as Kaiag 



Al VK^ (s-i) 6.68(13) 


5.83(11) 


1.20(12) 


9.82(10) 


5.19(7) 


VKg^ (s-i) 6.69(13) 


5.85(11) 


1.20(12) 


9.81(10) 


5.19(7) 


Sc VK^ (s-i) 3.96(14) 


2.50(13) 


7.02(13) 


3.96(11) 


1.39(9) 


l^si(s-i) 3.96(14) 


2.50(13) 


7.02(13) 


3.95(11) 


1.39(9) 



Table 3. Auger decay rates (ADR) per initial level (LSJi) of the Ti ground 
configuration 2s^ 2p^ 3s^ 3p^ 3d^ 4s^ obtained by adding the values of radiationless 
transition probabilities for all levels of the final configuration Is 2p® 3s^ 3p^ 3d^ 4s^. 



LSJ^ 


ADR (s-i) 


^Po 


1.592(14) 


^So 


1.592(14) 


^Pi 


1.594(14) 


^F2 


1.594(14) 


1D2 


1.594(14) 


^P2 


1.593(14) 


^F3 


1.594(14) 


^F4 


1.594(14) 


1G4 


1.593(14) 
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Table 4. Ratio of the K/3J' to K/3^ hypersatcUitc lines intensities, K/Si /Kf3^^ , computed 
in this work, with and without the Brcit interaction included, and compared with Chen 
et al 10 . 







K/3h/K/?i^ 








This work 




z 


Without Breit 


With Breit 


Chen 


18 


0.020 




0.017 


0.0093 


20 


0.040 




0.035 


0.022 


29 


0.32 




0.38 




30 


0.44 




0.45 


0.389 


36 


0.84 




0.88 


0.831 


38 


0.96 




1.00 




40 








1.10 


45 








1.34 


46 


1.33 




1.39 




47 








1.42 


48 


1.39 




1.46 




49 








1.48 


54 


1.55 




1.62 


1.61 


56 


1.59 




1.66 




60 








1.72 


65 








1.75 


70 


1.77 




1.82 




80 


1.84 




1.87 





Table 5. Values of widths for the Is ^ levels of Al and Sc calculated in this work. 





LSJ, 


r. (eV) 


Al 


^Pi/2 


1.079 




'P3/2 


1.022 


Sc 


'D3/2 


1.645 




'D5/2 


1.669 
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Table 6. K/3i and K/S^energy shifts, in eV, computed in this work and compared with 
Chen etal ^. 





S(K/?h)-^(K/3i) 


E{Kp^)-E{K(33) 


E{Kp^)-EiK(3^) 


z 


This work 


Chen 


This work 


Chen 


This work 


18 


226.04 


226.3 


225.34 


225.4 


0.88 


20 


255.22 


255.5 


254.37 


254.3 


1.24 


25 




324.8 




324.5 




29 


381.62 




381.14 




2.96 


30 


396.03 


397.4 


396.52 


397.0 


3.43 


31 


412.10 




412.15 




3.65 


36 


491.47 


492.3 


491.22 


492.1 


8.25 


38 


524.43 




524.17 




10.87 


40 




558.8 




558.6 




45 




645.7 




645.2 




46 


662.45 




661.86 




28.67 


47 




681.8 




681.2 




48 


698.88 




698.11 




35.50 


49 


717.42 


718.8 


716.56 


718.0 


39.22 


54 


813.56 


815.5 


812.14 


814.0 


63.60 


56 


854.28 




852.22 




76.43 


60 




941.7 




938.7 




65 




1055.6 




1051.1 




70 


1175.10 




1168.26 




229.67 


80 


1451.49 




1438.14 




444.43 
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Table 7. The K^^ ^ energy shifts calculated in this work and the available 
experimental results. 







K/3h3-K/3i,3 


energy shifts (eV) 


z 


This work 


Experiment 


13 




157.82 




14 




171.34 




15 




184.40 




16 




197.63 




17 




211.39 




18 




225.24 




20 




254.16 




21 




267.66 




23 




295.36 




25 






345 ± 35=^ 


26 




337.75 


336.0 ±0.5'^ 


27 




351.85 




28 




366.46 




29 




380.36 




30 




395.58 




31 




411.01 


390 ± 20== 


36 




489.78 




38 




522.60 




46 




660.00 




48 




696.26 




49 




714.25 


830 ± 60^= 


54 




810.42 




56 




850.82 




70 




1169.63 




80 




1443.46 




^ Ref. 


m 






''Ref. 


m 






" Ref. 


m 
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Table 8. Brcit and QED contributions to K(3i energy shift. 







Brcit 




QED 


z 


(eV) 


Percentage 


(eV) 


Percentage 


18 


1.66 


0.7% 


-0.08 


-0.04% 


20 


2.31 


0.9% 


-0.10 


-0.04% 


29 


8.59 


2.0% 


-0.30 


-0.08% 


30 


8.30 


2.1% 


-0.29 


-0.07% 


36 


14.77 


3.0% 


-0.47 


-0.09% 


38 


17.50 


3.3% 


-0.54 


-0.10% 


46 


31.96 


4.8% 


-0.88 


-0.13% 


48 


36.54 


5.2% 


-0.99 


-0.14% 


54 


52.79 


6.5% 


-1.61 


-0.20% 


56 


59.51 


7.0% 


-1.53 


-0.18% 


70 


121.72 


10.4% 


-2.96 


-0.25% 


80 


188.49 


13.0% 


-4.59 


-0.32% 
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Figure 1. Calculated spectrum of Si K/^Jg lines. A lorentzian with r=0.8 eV was 
used for each of the 14 transitions, which yield this profile. 
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Figure 2. Calculated K/Sj'g energy shifts compared with available experimental 
results. 



Relativistic calculation of Kj3 hypersatellite transitions 



17 



10%-, 



5%- 






0% 



-5% 



-10»/o 



-15% 



-20% 



-25% -I 



— 1 — 
10 



— I — 
20 



—I — 
30 



— I — 
40 



—I — 
50 



— I — 
60 



—I — 
70 



80 



Figure 3. Relative contribution of Breit interaction on the Kj3i/Kf3^ intensity ratio. 
BI and BN stand for Breit included and Breit neglected, respectively. 



